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Abstract 

Repulsive Bose-Einstein condensates immersed into a double-well trap potential are studied 
within the framework of the recently introduced mean-field approach which allows for bosons to 
reside in several different orthonormal orbitals. In the case of a one-orbital mean-field theory 
(Gross-Pitaevskii) the ground state of the system reveals a bifurcation scenario at some critical 
values of the interparticle interaction and/or the number of particles. At about the same values of 
the parameters the two-orbital mean-field predicts that the system becomes two-fold fragmented. 
By applying the three-orbital mean field we verify numerically that for the double-well external 
potential studied here the overall best mean-field is achieved with two orbitals. The variational 
principle minimizes the energy at a vanishing population of the third orbital. To discuss the energies 
needed to remove a boson from and the energies gained by adding a boson to the condensate, 
we introduce boson ionization potentials and boson affinities and relate them to the chemical 
potentials. The impact of the finite number of bosons in the condensate on these quantities 
is analyzed. We recall that within the framework of the multi-orbital mean- field theory each 
fragment is characterized by its own chemical potential. Finally, the stability of fragmented states 
is discussed in terms of the boson transfer energy which is the energy needed to transfer a boson 
from one fragment to another. 

PACS numbers: 03.75.Hh,03.65.Ge,03.75.Nt 
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I. INTRODUCTION 



The experimental realization of Bose-Einstein condensation provides an unique opportu- 
nity to study macroscopic quantum systems. One of the key property of a quantum system 
is the formation of shell structures. For a system of identical fermions the formation of 
the shell structure is enforced by the Fermi-statistic, which implies a restriction to the num- 
ber of particles residing in a single quantum level. In contrast, for bosonic systems where 
the Bose-statistic permits any occupation numbers, the formation of shell structures, also 
known as the fragmentation phenomenon must be caused by other reasons. 

Fragmentation can appear naturally in condensates made of different kinds of bosons. 
Binary mixtures of trapped Bose-Einstein condensates have been the subject of numerous 
experimental and theoretical investigations 14, nl l7l These studies comprise 

binary mixtures made of two different alkalis such as ^^i?6 — ol, of two different 
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same alkali such as the (F = 2,Mf = 2) and (F = 1,Mf = -1) states of ^'^Rb j^. 

In contrast, the fragmentation phenomenon in BEG made of atoms of one kind and in 
the same internal state, i.e., made of an identical bosons, is an open theoretical frontier. 
Traditionally, BEG is described at the Gross-Pitaevskii (GP) mean-field level ^, 1^, where 
all bosons are residing in a single one-particle state. This one-orbital mean-field has been a 
very successful approximation and can explain many experiments; see, e.g., Refs. jlll ll^ 
and references therein. However, this mean-field intrinsically cannot describe fragmentation. 
The best mean-field approach allowing for bosons to reside in several different orthonormal 

nn 

orbitals has been formulated recently [Ij, llJI . Each of the involved orbitals is character- 
ized by a different spatial distribution (localization in space), a different occupation number 
(number of bosons residing in the orbital), and a different chemical potential. This multi- 
orbital approach includes the Gross-Pitaevskii theory as a special case, namely when all 
bosons reside in a single orbital. The multi-orbital mean-field approach is based on a varia- 
tional principle and hence the optimal orbitals and their occupation numbers are determined 
variationally such that the energy of the system is minimized. In practice this means that 
the question whether a BEG forms a fragmented state (shell structure) or prefers to stay in 
a non-fragmented state (GP) is answered in the framework of one and the same method. 
By applying this multi-orbital mean-field to study a repulsive BEG in multi-well ex- 
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ternal potentials we found [15] that the ground state may be many-fold fragmented, i.e., 
the macroscopic occupation of several one-particle functions is energetically more favorable 
than the accumulation of all bosons in a single orbital. The quantities determining the 
fragmentation are the number of particles, the strength of the interparticle interaction and, 
of course, the specific shape of the external potential. The influence of all these parameters 
on fragmentation has been investigated in some detail. 

In this paper we consider a system of identical bosons with positive scattering length 
immersed into a double-well external potential. The Hamiltonian of this system is defined 
in Sec. II. In Sec. Ill we demonstrate that within the standard one-orbital mean- field theory 
(Gross-Pitaevskii), the ground state of the system reveals a bifurcation scenario starting 
from a critical value of the interparticle interaction strength or of the number of particles. 
By applying two-orbital mean-field theory we show in Sec. IV that the ground state of the 
system made of N identical bosons can indeed be two-fold fragmented. The fragmentation 
phenomenon starts to take place in the vicinity of the bifurcation point obtained at the 
GP mean-field level. In Sec.V we apply three-orbital mean-field theory and demonstrate 
that for the external potential studied here the overall best mean-field is achieved for two 
orbitals, i.e., the inclusion of a third orbital in the calculation does not improve the mean- 
field description and the lowest energy is obtained if only two orbitals are occupied. 

Sec. VI opens a second part of the present investigation, where we introduce characteristic 
quantities which allow us to distinguish, at least in principle, between fragmented and non- 
fragmented states of the condensates. In particular, by comparing the systems of and 

± 1 bosons we adapt quantities relevant to fermionic systems and define vertical and 
adiabatic ionization potentials and affinities for bosonic systems. Their relevance to an 
experimental observation of fragmented states is also addressed. In Sec. VII we introduce a 
boson transfer energy as the energy needed to transfer a boson from one orbital to another 
one and use it to analyze the stability of fragmented states. The discussion of the large N- 
limit of the boson ionization potential, boson affinity and boson transfer energy is presented 
in Sec. VIII. Finally, Sec. IX summarizes our results and conclusions. 
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II. THE SYSTEM AND HAMILTONIAN 



Our general intention is to consider a system of identical bosons interacting via a 
5-function contact potential W{fi — r*,) = Xo6{fi — fj), where is the position of the i-th 
)oson and the nonlinear parameter Aq is related to the s-wave scattering length of the bosons 



12| . The Hamiltonian of this system takes on the standard form 



N r 



1=1 



N 
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We denote as h{r) = T + V{f) the unperturbed one-particle Hamiltonian consisting of the 
kinetic operator T and the external potential V{r). 

In this work we specifically study bosons trapped in the one- dimensional double-well 
external potential shown in the insets of FigO Effectively the trap is obtained as an "inner" 
potential 

and an "outer" trap Vouter which consists of an infinite wall at x = 9.57r where the inner 
potential has already died off. The resulting combined potential {V{x) = Vinner + Vouter) has 
two well-separated nonequivalent wells. It should be mentioned that the results discussed 
here hardly de pen d on whether the infinite wall of Vouter is replaced by a smoothly growing 
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potential wall [15||- The kinetic energy reads T = — implying that the coordinate x is 
dimensionless while all energies and Aq are now in units of the frequency u. 

We would like to stress that all conclusions and qualitative results discussed in this 
work also apply to other double-well external potentials like, for instance, that used in the 
experimental set up of Ref.jl^. 



III. GROSS-PITAEVSKII MEAN-FIELD RESULTS 

The standard one-orbital mean-field description of the interacting system is obtained by 
assuming that the ground state wave function \1/ is a product of identical spatial orbitals 
ip: \l/(ri, r2, . . . , -nv) = V^(^i)v^(^2) ■ ■ ■ V'(^Ar)- The energy E =< \I^|f/'|\I' >, defined as the 
expectation value of the H, reads 

Egp = N{J ^*h^dr + ^ j \v\'dr}, (3) 
4 



where A = Ao(A^ — 1) is the interaction parameter. By minimizing this energy the well-known 



Gross-Pitaevskii equation j9 



10| is obtained 



{ h{r) + Ao(iV - l)|^(rl|2} ^(r) = ficp^if^- 



(4) 



The only parameter involved in the Gross-Pitaevskii mean- field is A = Ao(A^ ^1)- There- 
fore, all systems which are characterized by the same A have the same energy per particle as 
well as the same orbital V9(r) even if the systems have different numbers of bosons. Below, 
we will use this fact to compare systems made of different numbers of bosons but have the 
same value of A. 

We solved the Gross-Pitaevskii equation with the external potential introduced in Sec. II 
for different values of A. In Fig^we plot the energy per particle of the lowest energy solutions 
as a function of A. From this figure it can be inferred that the energy per particle of the 
ground state increases monotonically with A up to some critical value of Xcr = 0.837 and then 
the energy trajectory is split into two branches. This critical value of A is indicated in Fig^ 
by a vertical dashed line. The single-particle wavefunction (orbital) corresponding to the 
upper branch, marked as I in Fig^ is depicted in the upper inset of FigQ This wavefunction 
is mainly localized in the left well. From Acr on, the lowest energy branch, labeled II in Fig^ 
smoothly bifurcates from the " localized" solution and describes a delocalization of the BEG 
over the two wells. This delocalized solution is depicted in the right bottom inset of FigO 
For convenience, we also plot in both insets the rescaled double- well potential 1^/10. The 
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of the Gross-Pitaevskii equation 



bifurcation scenario is a known feature 
refiecting the non-linearity of the un derlying mean-field approximation. Bifurcations have 



)een predicted for the ground 
3, Q 20, 21 1 and attractive 
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li La y, y, 



and excited [1^, |21| states of repulsive 



22| | condensates confined in different trap 



potentials. 

Within the Gross-Pitaevskii mean field we obtain the following physical picture of Bose- 
Einstein condensation of repulsive bosons immersed into non-symmetric double-well poten- 
tials. For A < Acr the BEG is localized in the deeper well only. When A is increased, the BEG 
may continue to be localized in this well, but it is energetically more favorable for bosons 
to tunnel through the barrier and to populate also the other well. As we have mentioned 
above, A = Ao(A^ — 1) is the only relevant parameter, and therefore both localization and 
delocalization phenomena can be observed for systems made of any number of particles by 
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properly tuning the interparticle interaction strength Aq. 



IV. TWO-ORBITAL MEAN-FIELD RESULTS 

In the previous section we demonstrated that above some critical value of nonlinearity 
Acr the ground state of the system is described in the GP approach by a wave function 
delocalized over both wells. The existence of a bifurcation in the GP results can be seen as 
a hint to go beyond the GP theory. The natural question arises whether it is possible to 
improve the GP mean-field description by providing a more flexible mean-field ansatz which 
would allow bosons to occupy two different orbitals which, in principle, could be localized 
in different wells. 

Recently, a multi-orbital mean-field approach allowing for bosons to reside in several 
different orthonormal one-particle functions has been formulated jl^l . Il^ . Using m different 
orbitals, the approach is denoted MF(m). In this section we discuss MF(2). In the next 
section we shall see that for the double-well geometry of the external potential studied here, 
the overall best mean-field is achieved within the two-orbital mean-field theory. Since the 
details of the derivations have already been published elsewhere jisl. [l^. we present here 
only the final formulae. Assuming that rii bosons occupy the orbital (f)i and ^2 bosons 
occupy the orbital 02 leads to the following many-body wave function 

*(fi, . . . , f^) = 50i(fi) ■ ■ ■ (f)i{fnj(f)2{rm+i) ■ ■ ■ 02(r„i+„2), (5) 

where S is the symmetrizing operator. The MF(2) energy expression takes on the form 



+2Aoriin2 j \(t)i\^\(t)2?dr (6) 



The optimal orbitals 0i and 02 which minimize this energy functional are determined by 
solving the two coupled non-linear equations: 

{ h{r) + Ao(ni - 1)101 (r-)|2 + 2Aon2|02(f)p} 0i(f) = 

= /Wll0l(^ +/il2 02(^ 

{ h{r) + AoK - l)|02(r-)|2 + 2Aoni|0i(f)p} 02(f) = 

= Ai2i 01 (r) + /^22 02 (r) (7) 



Obviously, the GP equations © and (jH) follow immediately from the MF(2) equations 
(0) and ((7j) by putting the occupation of one of the two orbitals to zero. 

In contrast to the one-orbital mean-field (GP) description where A = Ao(A^ — 1) is the only 
parameter involved, the two-orbital mean-field (MF(2)) depends on two parameters rii and 
n2 which are the occupation numbers of the one-particle orbitals 0i and 02 respectively. At 
fixed interaction strength Aq and number of bosons A^, the MF energy (see EqlHl) depends on 
the particular value of ni {n2 = N — rii). It should be noted that the occupation numbers are 
variational parameters. In order to find their optimal value and that of the energy, Eqs.((7)) 
are solved for all possible occupation numbers, and the corresponding energies are evaluated 
by using Eq.(jni). At the optimal values of the occupation numbers the energy E takes on 
its minimum and the MF(2) becomes the best two-orbital mean field which we denote by 
BMF(2). 

We solved the MF(2) equations for A = 0.8 and for A = 0.9. These values of the non- 
linear parameter have been chosen to be smaller and larger than Acr = 0.837 which is the 
bifurcation point of the GP energy trajectory discussed in Sec. III. In FigEland FigOwe plot 
the MF(2) energies per particle for A = 0.8 and for A = 0.9, respectively, as a function of 
the relative occupation number rti/N oi the orbital localized in the deeper (left) well of the 
trap potential. The results are shown for A^ = 5, 10, 25, 10^, 10^ and 10^ bosons. The GP 
energy is indicated by the horizontal line and labeled as " GP" . 

The common feature seen in Figs|2l and El is that the energies per particle of the systems 
of bosons which are characterized by the same A = Ao(A^ — 1) do depend on the number 
of bosons A^ and reveal a different dependence on rii/N . In contrast to the GP approach, 
the multi-orbital mean-field theory distinguishes between systems with different number of 
bosons even if they are characterized by the same value of A. As seen in the figures, for 
rii/N = 1, i.e., when all bosons reside in one orbital, the energies of all the curves coincide 
at the same value which is nothing but the corresponding GP energy. We can obviously 
conclude that the GP approach is indeed a special case of the MF(2) theory. 

In FigI21and in particular in FigOlit is seen that the optimal BMF(2) energies may be 
lower than the corresponding GP one and in these cases BMF(2) constitutes the proper 
ground state of the system. These ground states with nonzero occupation numbers are 
fragmented states, because the respective one-particle density matrix has several nonzero 
eigenvalues Q, |2^ . Therefore, the fragmentation observed before in three- well potentials 
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[15| also persists in double-well external potentials. It is worthwhile to stress that at the 
minimum of the energy one of the corresponding optimal orbitals is localized in the left well 
and the other one in the right well of the trap potential. This is in an agreement with the 
three-fold fragmentation studied before, where the optimal orbitals have also been found to 
be localized in the different wells of the three-well external potential |l5| . 




Figs|21 and El exhibit several interesting features. Let us first discuss FigO where A = 0.8 
is smaller than Acr- Inspection of the inset shows that the MP (2) energy curves as a function 
of rii/N posses a minimum for < 14 bosons. At these minima the energy is below the 
GP energy which is at the maxima of the curves. Por the systems which contain more than 
14 particles, the MP(2) energy curves are monotonously increasing and their lowest energy 
is at the GP energy. 

The situation changes drastically for A = 0.9 which is larger than Acr- As seen in the 
inset of PiglSl the BMP(2) energy is lower than the GP energy for condensates with up 
to Nmax ~ 102000 bosons. Por all these systems the fragmentation of the ground state is 
evident and macroscopic. One should be aware that A is kept fixed in PigOland, therefore, 
the interparticle interaction strength Aq is small for the many-particle systems shown there. 

Pinally, we would like to point out two additional features of Pig 01 Pirst, the GP energy 
is a local minimum for most of the energy curves and consequently there appears a maximum 
which separates this energy from the global minimum, i.e., from the BMP(2) energy. This 
may hint to an interesting dynamics of fragmentation if the system is initially prepared in an 
non- fragmented state. Second, although there exists a maximal number of particles, Nmax, 
for which the ground state is fragmented, it should be stressed that for all N > Nmax the 
energy gap between the BMP(2) and GP energies is extremely small. Por = 10^ bosons, 
for instance, this energy difference per particle is only ^ 3 x 10~^ units of u. The coexistence 
of the fragmented and non-fragmented, but delocalized states is of great interest by itself 
and might play a role in particular in time-dependent experiments. 

V. PROOF THAT BMF(2) IS THE BEST OVERALL MEAN FIELD 

Prom the previous section we learned that by allowing bosons to reside in two orbitals, 
the mean-field description of the repulsive BEG in the double-well potential has improved 
compared to the standard one-orbital description. The natural question arises whether the 
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inclusion of even more orbitals can provide further improvement of the mean-field description 
of the BEC. To answer this question we apply the three-orbital mean-field theory (MF(3)) 
[l5|. 

The corresponding ansatz for the wavefunction assumes that three orbitals 0i, 02 and 03 
are now occupied by ni,n2 and = N — rii ~ n2 bosons, respectively 

^{n, . . . , rW) = 50i(fi) • ■ ■ Mrm)Mrm+l) ■ ■ ■ (j)2{rm+n2)Mrm+n2+l) " " " 03(n^l+n2+r^3), (8) 



where S is the symmetrizing operator. The three-orbital mean-field energy reads 



+2Aonin3 j 101^031^^?^ + +2 Aongns j |02n03pc?^- (9) 

The optimal orbitals minimizing this energy functional are obtained by solving the following 
system of three coupled non-eigenvalue equations [3| 

{ h{r) + Xo{n, - 1)101 (f)|2 + 2Aon2|02(r)|' + 2Aon3|03(r-)p} 0i(f) = 

= yUll 01 (r) + flu 02 (r) + fll3 03 (r) 

{ h{r) + Ao(n2 - l)|02(r1p + 2Aoni|0i(f)|2 + 2Aon3|03(r-)p} 02(r-') = 

= ^21 (Pl{r) + /i22 02 (^) + yU23 03(0 

{ h{r) + AoK - l)|03(r)P + 2Aoni|0i(f)|2 + 2Aon2|02(r-)p} 03(r) = 

= yU31 01 (0 +/^32 02(0 +yU33 03(0- (lO) 

For the numerical procedure to solve this system of equations and to obtain the self-consistent 
orbitals 0j and the corresponding values of the Lagrange parameters /ijj = 1,2,3), we 
refer to Ref . jisf . 

In comparison with the MF(2) approach the MF(3) method has one more variational 
parameter n^. We solved the above system of three coupled equations for several values 
of A and and different fixed values of n^. As an example, we show in Fig^the results 
obtained for A = 0.9 and = 1000 for different occupation patterns. Keeping the value of 
the third occupation number n^/N fixed at 0.0001, 0.001, 0.005 and also at zero, we plot the 
three-orbital mean-field energy per particle as a function of the relative occupation number 



9 



rii/N . The occupation of the second orbital n2 = N — ni — is, of course, determined 
by 111 and n^. The three-orbital mean- field method obviously includes MF(2) as a special 
case, namely, when the occupation of the third orbital vanishes. We see from FigjU that 
increasing the value of n^/N leads to a gradual increase of the MF(3) energy per particle. 
Clearly, in the present case of repulsive condensates in the double-well potential, the overall 
best mean-field is obtained within the two-orbital BMF(2) theory. The enforced inclusion 
of more orbitals only enhances the energy of the condensate. 

Several consequences should be mentioned. First, fragmentation is a general physical 
phenomenon which takes place in repulsive BECs trapped in multi-well external potentials. 
Second, the number of the fragments and their occupation numbers are defined variationally, 
by minimizing the total energy functionals in Eqs.® and Q. If more orbitals are included 
in the mean-field ansatz than needed, the occupation of the superfluous orbitals becomes 
zero. 



VI. SOME PROPERTIES OF FRAGMENTED STATES 



Usually, to distinguish fragmented and non-fragmented states a phase difference between 
wave functions of the different fragments is considered [2^, U^- Dynamical stability 



271 128| of the relative phase between the condensates localized in the different wells and 
related questions on the evolution of the fragmented state 3, 2^ have been a subject of 
several discussions. 

In contrast to these time-dependent studies on fragmentation, we concentrate in the sec- 
ond part of our work on time-independent properties of the stationary fragmented states 
themselves. The analogy between shell structures (atoms and molecules) formed in the 
fermionic world and the bosonic fragmented states studied in the previous sections moti- 
vates us to make use of some fermionic observables such as the ionization potential and 
electron affinity, and introduce related quantities for bosonic systems. Indeed, at the multi- 
orbital mean-field level of description the fermionic system is described by antisimmetrized 
single-determinant wave function [s^, while the symmetrized ansatz (see Eqs.© and (jH))) 
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is proposed for bosonic systems. The variationally optimal orbitals and the correspond, 
orbital energies are obtained by solving the well-known Hartree-Fock (HF) equations 
for fermions and the BMF equations for bosons. The total energies of these systems are 
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evaluated by computing the expectation values of the full many-body Hamiltonians with 
the (anti) symmetrized ansatz of the wave function. 

The electron ionization potential (IP) is defined as the difference in total energies between 
the reference system with electrons and the ionized system with — 1 electrons jsi]. 
Because of this definition, the ionization potential is also referred to as the binding energy 
of an electron in the system. We define the boson ionization potential or the binding energy 
of a boson as the energy needed to remove a boson from a bosonic system with bosons. 
Another very important physical characteristic is the electron affinity defined as the energy 
gained by adding one electron to the atom or molecule 0|. In analogy, we define the 
boson affinity as the energy gained by attaching one boson to the bosonic system under 
consideration. 

The boson ionization potentials and affinities can be calculated straightforwardly based on 
their definitions as differences between total energies. In principle, these total energies may 
be evaluated within the framework of any suitable N-body method, but in the present study 
we restrict ourselves for consistency to the two-orbital mean-field approach. To distinguish 
between the total energies of the condensates made of a different number of bosons we 
introduce superscript indices. In this notation the MF(2) total energies of the system with 
A^— 1, A^ and A^-|-l bosons take on the following form '^~^E{ni,n2), '^E{ni,n2), '^^^E{ni,n2), 
where the superscript refers to the condensates made of the corresponding number of bosons. 
We recall that within the framework of the MF(2) approach the total energy of the system 
is a function of the occupation number ui {n2 = N — ui) and the minimum of this energy 
is called the best mean field (BMF(2)). Obviously, the MF(2) energies of the systems with 
a different number of bosons achieve their minimal values (BMF(2) energies) at different 
values of the occupation numbers. We denote these BMF(2) energies and corresponding 
occupation numbers as ^Eq and {nf , nf^) where the superscript K refers to the system with 
K bosons. 

A. Adiabatic and Vertical boson ionization potentials 

The general definition of the boson ionization potential involves the total energies of 
condensates made of A^ and A^ — 1 bosons. Usually, the reference system with A^ bosons is 
considered to be in its ground equilibrium state while the ionized system with A^ — 1 bosons 
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can be either in an equilibrium or in a some transitional state. In the following section we 
verify that only the BMF(2) energy (the minimum of the MF(2) energy curve) corresponds 
to the equilibrium state of the condensate, while all other MF(2) energies can be attributed 
to some transitional (non-equilibrium) states of the condensate. 

The difference between the ground equilibrium state energies of the ionized system with 
— 1 bosons and the reference system with bosons is called adiabatic boson ionization 
potential (BIPa)'- 

BIPa = - ^i?o. (11) 

The equilibrium in the ionized state might be achieved by adiabatically removing a boson 
from the condensate. 

In contrast, if a boson is suddenly removed from the system, then the state created is 
not an equilibrium state of the ionized system and the difference between the energy of this 
non-equilibrium state and that of the ground state of the reference system is called vertical 
boson ionization potential (BlPy). If the bosonic system under consideration is fragmented, 
then the sudden ionization of a boson from different orbitals of the fragmented state requires 
different energies. The two-fold fragmented states are characterized by two different vertical 
boson ionization potentials 

BlPvil) = - - ^i?o«,n^) 

BIPy{2) = ''-'E{n^,n^ - 1) - ^i?o«,0, 

(12) 

where '^Eq = '^EqIu'^ ,n2) is the BMF(2) energy of the reference system with bosons 
obtained at the optimal occupation numbers and of the two fragments and, ^~^E{n^ — 
1, n^) and ^~^E{ni , — 1) are the MF(2) energies of the ionized states where a boson has 
been removed from the first and second fragment, respectively. 

In FigElA. we present a schematic diagram of the adiabatic and vertical ionization poten- 
tials together with the MF(2) total energies of the systems with + 1, and A^ — 1 bosons 
plotted as functions of the number ni of bosons residing in the first orbital. In Fig|33 we 
plot these quantities as functions of the complementary parameter n2 = N — Ui (the number 
of bosons residing in the second orbital). The BMF(2) energies of the systems with A^ + 1, A^ 
and A^ — 1 bosons are indicated on the energy axis. The optimal occupation numbers 



12 



and corresponding to the minimum of the MF(2) energy of the reference N-boson system 
are indicated on the x-axes of Figl^A. and Fig|33, respectively. 

The difference between the ground state energy of the ionized system ^~^Eq and that 
of the reference system ^Eq, i.e., the adiabatic ionization potential, is marked in FigEK as 
BIPa. The vertical boson ionization potential BIPy{2) is shown in FigEJA. as a vertical 
line connecting the minimum of the A^-particle energy curve and the energy curve of the 
system with — 1 bosons. Let us explain this construction. When a boson is suddenly 
removed from the second orbital 02 of the reference system, the occupation number of 
the other (first) orbital remains the same also in the ionized system, while the occupation 
of the second orbital in the ionized state is obviously reduced by one to — 1. Therefore, 
this ionization process is indicated by the vertical line in the (E', ni)-diagram. 

The vertical boson ionization potential BIPy{l) is also shown as a vertical line connecting 
the minimum of the N-particle energy curve and the energy curve of the system with — 1 
bosons, but in the complementary {E,n2 = N — 'rii)-diagram (see FiglH^). In this case the 
occupation number of the first orbital is reduced by one to — 1 in the ionized state, while 
the occupation number of the second orbital remains the same as in the reference state. 
For completeness, the BIPy{l) and BIPy{2) are plotted in both panels of figure 

Removing a single boson from the system will not cause a strong change of the orbitals 
of the other bosons and these orbitals can be assumed fixed. In fermionic systems this 
assumption is known as the frozen orbital approximation [31]. This observation allows us to 
evaluate approximately the vertical ionization potentials and hence, the total energies of the 
ionized systems. 

In the frozen orbital approximation we assume that the orbitals of the reference N-bosonic 
system 0i and 02 do not change upon the sudden removal of a boson, i.e., the ionized system 
is described by the same orbitals as well. The energies of the corresponding ionized states 
where the boson has been suddenly removed from the 0i or the 02 orbital take on the form 
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^-iE(ni-l,n2) = (m - l)/^ii + Ao^^^^^l^^^ I \M'df+ 

n2/i22 + Ao^^^^^^l^ j |02|'rfr + 2Ao(ni- 1)722 j \<pi\^\<p2\^df (13a) 

N-lr.,^ „ 1X / , X (n2-l)(n2-2) /• 4 



-E(ni,n2-1) = (r22 - l)/i22 + Ao^^ J |02|W+ 

ni/iii + Ao^^^i^^^P^ j |0i|^c/f++2Aoni(n2-l) ^ |0iri02rc?r (13b) 

respectively. 

According to the definitions in Eqsll2| the vertical ionization potentials are obtained by 
substructing the BMF energy '^Eq of the reference N-particle state from the energies of the 
ionized states in Eqs ]13all3bl at n2 = n^. We easily find for the boson vertical ionization 
potentials in the frozen orbital approximation the appealing relations 

BIPvil) = -/in 

BIPvi2) = -/i22, (14) 

i.e., they are given by the negative of the respective Lagrange multipliers which we call 
chemical potentials. 

It is worthwhile to recall that at the level of the standard GP theory (see Sec lIII|) the 
energy needed to remove a boson from the non-fragmented N-boson condensate without 
changing the corresponding orbital ip is given by the GP chemical potential: 

^~^Egp — '^Egp = —fJ-GP- (15) 

The relations between orbital energies and ionizati on p otentials in the fermionic case 
are the subject of the well-known Koopmans' theorem |3ll. [3^. Whithin the Hartree-Fock 
mean field this theorem states that the vertical ionization potential is the negative of the 
energy of the orbital from which the electron has been removed. Attributing the diagonal 
Lagrange multipliers (chemical potentials) fin to the orbital energies of the fragmented state, 
the formal results for fermionic and bosonic systems are absolutely identical. 

We conclude that the process of sudden ionization of a boson from the non-fragmented 
GP state is characterized by a single ionization potential while different ionization potentials 
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characterize a two-fold fragmented state of the condensate. We recall that in the thermody- 
namic limit the change in the energy of a system obtained by removing a particle is called 
chemical potential. The condensates studied here are made of identical bosons, and, there- 
fore, a single chemical potential is expected to characterize these condensates at equilibrium. 
This seems to contradict the results obtained here for the two-fold fragmented states, where 
two, in general different, chemical potentials exist. This aparent contradiction is resolved 
by remembering that the present results are for finite systems. Indeed, we demonstrate in 
Sec JVIIll that in the limit of a large number of particles these different chemical potentials 
become identical at the optimal occupations restoring thereby the thermodynamic picture 
of a condensate. 



B. Adiabatic and vertical boson affinities. 



The main propose of the present section is to consider a process where a boson is added 
to the reference system of bosons. We define adiabatic boson affinity {BAa) as the 
difference between the ground equilibrium state energies of this reference system and of the 
system which results by attaching a boson to the reference system 

BAa = ^i?o - (16) 

The adiabatic boson affinity is shown schematically in Fig^A.. It describes the energy gained 
by the attachment of a boson. 

Similarly to the sudden ionization process used to introduce the vertical ionization po- 
tentials in the previous subsection, we consider the sudden attachment of a boson to the 
two-fold fragmented state and define two vertical boson affinities comprising the sudden 
attachment of a boson to the 0i and 02 fragments, respectively: 

BAyi2) = ^Eo(nf,0- + 

(17) 

Here, ^Eq = -^£"0(71^,71^) is the BMF(2) energy of the reference system, and are 
the corresponding optimal occupation numbers, ^~^^E{n^ + ^,1^2) and ^^^E{n^ + 1) 
are the MF(2) energies of the states where a boson has been attached to the first and to the 
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second fragment, respectively. In Figl^A. and Figl^fB we schematically plot the vertical boson 
affinities BAv{2) and BAv{l) as the vertical lines connecting the corresponding points on 
the energy curves. 

In the frozen orbital approximation the energies of the states to which the boson has been 
suddenly attached to the 0i or 02 orbital take on the form 



^2/^22 + Ao^^^%^ / |02|'cir + +2Ao(ni + 1)712 / Ihl^lhl^ df (18a) 



^+^i?(ni,n2 + l) = (n2 + l)/.22 + Ao^^^^^^l^^ J \h\'df+ 
r^i^ii + Ao^^^i%^ / |0i|'rfr + +2Aoni(n2 + l) [ \<f)i\^\(f)2\^df (18b) 



2 

By substructing these energies from the BMF energy of the reference N-particle state ^Eq, we 
obtain the vertical boson affinities BAv{l) and -8^4^(2) in the frozen orbital approximation: 

BAv{l) = -/iii-Ao j \(pi\^df 

BAvi2) = -/X22-A0 J |02|'c^r. (19) 

Interestingly, the vertical boson affinities in the frozen orbital approximation are not fully 
determined by the chemical potentials in contrast to our finding for the vertical boson 
ionization potentials discussed in the preceding subsection. 

It is informative to notice that the energy gained by adding a boson to the non-fragmented 
condensate without changing the corresponding orbital if for a finite number of bosons reads: 

""Egp - ""^'Egp = -fiGP - Ao y \^\'df. (20) 

Hence, even at the level of the standard GP theory, the energy needed to remove a boson from 
a non-fragmented condensate differs from that gained by adding a boson to this condensate. 

The above findings reveal differences between bosonic and fermionic systems. In contrast 
to the electron affinity which in the framework of the Hartree-Fock mean-field approach is 
given by the negative of the virtual orbital's energy, the boson affinity does not depend on 
the virtual orbital and, in addition, is subject to a correction term proportional to Aq. 
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C. Numerical examples 



For illustration purposes we evaluate all the above introduced adiabatic and vertical 
boson ionization potentials and boson affinities for several reference systems made of = 
5, 10, 25, 100 and 1000 bosons, keeping A = 0.9 fixed throughout. In the GP theory all 
these systems are characterized by the same chemical potential. Tab]!] summarizes the 
computational results on the vertical boson ionization potentials and vertical boson affinities 
obtained within the framework of the direct scheme in Eqs^l and El respectively. For 
comparison we also evaluate all these quantities using the frozen orbital approximation as 
given in EqsElaiidEl In Tab|l] these approximate numbers are given in parenthesis. By 
comparing all these quantities we conclude that the frozen orbital approximation provides 
very accurate results for the systems with a large number of bosons and surprisingly accurate 
results for the system with a small number of bosons. 

The main physical conclusion derived is that the two-fold fragmented states are char- 
acterized by two different vertical ionization potentials and by two different vertical boson 
affinities. In contrast, the non- fragmented GP state is characterized by a single ionization 
potential and single affinity. At least in principle, this difference might be used to distinguish 
fragmented and non-fragmented states. The values of the various vertical boson ionization 
potentials and vertical boson affinities differ from each other clearly for the systems made 
of a small number of bosons. As the number of bosons increases, all these values become 
more similar to each other and also to the negative of the GP chemical potential. 

Closing this section, we would like to mention that we used the numerical MF(2) data 
obtained for the reference system with N = 25 and for the A^ ± 1 systems with 24 and 26 
bosons to plot the diagrams presented in FigsEl For all these systems we fixed the interaction 
strength at Aq = 0.9/(25 — 1). The optimal occupations of the ffist and second orbitals at the 
minimum of the MF(2) energy curve of the reference system with 25 bosons are found to be 
= 23.08 and ra^ = 1.92, respectively. The non-integer occupation numbers appear due 
to the underlying mean-field approximation. They are to be considered as average values of 
the respective particle-number operators. We use these occupation numbers as the reference 
scale for the x-axes in Fig|3lA and Fig|2|3. 
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VII. BOSON TRANSFER ENERGIES AND THE ORIGIN OF THE MINIMA OF 
MF ENERGIES 



We define the boson transfer energy (BTE) as the energy needed to transfer a boson from 
one orbital to the other. In the framework of the MF(2) theory we start with the system 
which has the energy ^E{ni, ^2) and evaluate the BTE as the energy difference between this 
energy and the energy of the "final" state where the occupation number of the first orbital 
is reduced by one to rii — 1 and the occupation number of the second orbital is increased by 
one to n2 + 1. In principle, the energy needed to move a boson from the first orbital to the 
second one differs from the energy of the inverse process where a boson is transferred from 
the second orbital to the first one. The BTE(i^2) and BTE(i^2) are given by: 

STE(i^2) = ^£^(^1 T 1,^2 ± 1) - '^^(ni,n2). (21) 

In the frozen orbital approximation the boson transfer energies can be straightforwardly 
evaluated using these definitions and the MF(2) expressions for the energies. The results 
take on the following form: 

5T^(i^2) = /i22-/^ii + Y j |02|'rfr-2Ao j \<pi\''\<p2?dr 
BTE^,^2) = /Xn-/X22 + y j \<pi\'df ~ 2 j \<pi\'\<p2\'dr 

(22) 

These results will be used in the discussion below. 

It is quite natural to suppose that at equilibrium an exchange of the particles between 
different orbitals (fragments) must be suppressed in the fragmented state. In other words, 
the energy needed to transfer a boson from the first orbital to the second one, is equal to 
the energy needed for the inverse process - a transition of the boson from the second orbital 
to the first one. Therefore, the difference between both boson transfer energies 

Q = (5TE(i^2) - BTE^,^2))/2 (23) 

may serve as a criterion of how far we are from the equilibrium. Substituting the respective 
transition energies from Eqs]22l provides a senseful approximation for Q: 

Q = /ill - ^^22 + y y 101^^^- ^ j \<P2\'df (24) 
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To make contact with the MF(2) energy curves and to get deeper insight into the origin of the 
minima of these curves which determine the BMF(2) energy, it is worthwhile to consider the 
energy derivative with respect to the occupation number rii. To evaluate the derivative 
we use again the frozen orbital approximation which has been found to provide accurate 
results for the boson ionization potentials and affinities. In this approximation we can 
assume that the orbitals 0j do not depend explicitly on the and the direct differentiation 
of the energy in EqlHlwith respect to rii gives: 

^ = //II - /i22 + y J |0i|^rfr-y J \(l)2\'^df (25) 

By comparing this derivative with the previous equation for Q one finds that they are 
identical. 

At the optimal orbitals and occupation numbers the MF(2) energy takes on its minimum 
and ^ = (see FigsO and Ej). Consequently, our physical assumption that at the equi- 
librium the boson transfer energies are equal is fully supported by identifying Q with 
Moreover, if the bosonic system is not at equilibrium, then due to the difference between 
BTE(i^2) and BTE(i^2) a fiow of bosons between the 0i-manifold and the 02-nianifold is 
enforced until equilibrium is reached. Therefore, the slope of the MF(2) energy curve given 
in Eq|2nican be viewed as the "driving force" for the fiow of bosons between the two boson 
subsystems. 

To verify these physically appealing results we investigate in more detail condensates 
with = 25 and = 1000 bosons at A = 0.9. The lower panels of FiglHlshow values of 
Q for different values of rii/N evaluated via Eq|21] using the frozen orbital approximation. 
Also shown is ^ obtained by numerically differentiating the MF(2) energy curve in FiglHl 
with respect to ni/N for the same systems. On both pictures the values of the optimal 
occupation numbers rii/N where the MF(2) energy curves take on their minima are marked 
by vertical dashed lines. From these figures it is clearly seen that indeed, at the minimum of 
the MF(2) energy, the fiow of the particles between the fragments is completely suppressed, 
i.e., at the best mean-field state the driving force is equal to zero. 
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VIII. LARGE N LIMIT 



In the large N limit, as well as ni and n2 are much larger than 1 and we may replace 
— 1 by A^. Repeating the calculations of the preceding sections one finds that all terms 

proportional to Aq x 1 and not to Aq x or Aq x nj, i = 1, 2, vanish in the large N limit. 

Accordingly, we find that in the large N limit not only the vertical boson ionization potentials 

are given by the chemical potentials as found in Sec lVII fsee Eqs J14|) . but that also the vertical 

boson affinities are: 

BAvik) = -fikk, k = l,2. (26) 

Obviously, also the boson affinity in the framework of the GP theory is determined by the 
corresponding GP chemical potential ficp in the large N limit. 

In the large N limit the boson transfer energy BTE(i^2) becomes identical to —BTE(^i^2)- 
The energy needed to transfer a boson in a condensate with A^ particles from one fragment 
to the other one is just given by the difference of the corresponding potentials: 

BTEii^j) = i^jj - jjii (27) 

This quantity is then nothing but BIPy{i) — BAy{j) which is very appealing for large 
systems. 

In the large N limit the quantity Q introduced in EqI2Sl and the energy derivative with 
respect to the occupation number ni become identical and are just determined by the dif- 
ference of the chemical potentials of the fragments. In particular, we find 

dE . . 

— = /ill - /^22 (28) 

aril 

which has been derived before in Ref . jl^ . 

To illustrate the impact of the growing number of bosons in condensates we plot the 
diagonal Lagrange multipliers /in and ^22 as a function of the relative occupation number 
rii/N for systems with A^ = 25 and A^ = 1000 bosons in the upper panels of FigEl In 
these panels the values of the optimal occupation numbers ni/N obtained at the minima 
of the MF(2) energy curves E{ni) are marked by vertical dashed lines and the value of the 
chemical potential ficp obtained in the GP theory by a horizontal solid line. 

One can see that for the systems with A^ = 25 bosons (see left upper panel of FiglHI) the 
crossing point of the fj,ii{ni/N) and fJ^22{ni/N) curves is located substantially aside from the 
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vertical dashed line. In contrast, the corresponding curves cross the vertical dashed line at 
almost the same value oi ni/N for the system with = 1000 bosons (see right upper panel 
of FigEI). However, on an enlarged scale as seen in the inset, the exact crossing point of the 
fJ'ii{i^i/N) and fi22{i^i/N) curves lies aside from the vertical dashed line marking the value 
of the optimal occupation number ni/N. 

To better understand these results let us analyze Eg 031 The systems with = 25 
and = 1000 bosons at A = 0.9 are characterized by different values of Aq = A/(A^ — 1). 
Therefore, the contributions of the terms which are proportional to Aq in Eql^ are much 
smaller for the system with = 1000 than with A^ = 25 bosons where they play an 
important role. In other words, in the large N limit the flow of the bosons is determined by 
the differences of the chemical potentials /in — /i22 only. Interestingly, = holds at the 
best mean field and consequently the chemical potentials of the fragments are identical to 
each other in the ground state of the condensate in the large N limit. 

IX. CONCLUSIONS 

In this article condensates with A^ identical repulsive bosons immersed into a double-well 
external potential have been investigated at different levels of multi-orbital mean-field the- 
ories where optimal orbitals and optimal occupation numbers are determined variationally. 

At the level of one-orbital mean-field theory (Gross-Pitaevskii) the ground state of the 
system reveals a bifurcation scenario at a critical value of the non-linear parameter Acr = 
Ao(A^ — 1) = 0.837. From this Acr on the delocalization of the bosons over the two wells of 
the trap potential becomes energetically more favorable than the localization of the BEG in 
the deeper well. 

The two-orbital mean-field theory predicts the existence of two-fold fragmented states, 
i.e., both orbitals of these states are occupied macroscopically. Depending upon the number 
of particles in the BEG and/or the strength of the interparticle interaction, this fragmented 
state can be the ground state of the system. For condensates with a small number of 
particles the fragmented ground state is the only stable state, because the non-fragmented 
state predicted by the GP theory corresponds to the maximum of the energy in the two- 
orbital scenario. 

By applying the three-orbital mean-field we verify numerically that for the double-well 



21 



external potential studied here the overall best mean field is achieved within two orbitals 
(BMF(2)). The inclusion of third orbital does not improve the mean-field description as 
the variational principle minimizes the energy at a vanishing population of the third orbital. 
Hence, the two-fold fragmented state obtained is the physical state of the system. Once N > 
Nrnaxi whcrc Nmax dcpcnds on the interparticlc interaction strength, the energy difference 
between the fragmented and the non-fragmented state may be very small, indicating that 
in these cases two physically different states may coexist. 

Inspired by similarities between fragmented states in bosonic systems and shell structures 
in the fermionic world, we introduce vertical and adiabatic boson ionization potentials and 
boson affinities. In this respect we discuss the chemical potentials of the fragments and their 
relation to boson ionization and attachment measurements. The energy splitting between 
the values corresponding to the different fragments are small for systems with a large number 
of bosons while they might be distinguishable for systems with a small number of bosons. 

We also introduced a boson transfer energy as the energy needed to move a boson from 
fragment " 1" to fragment "2" . We proved that only at equihbrium the same energy is needed 
for the inverse process where a boson is transferred from fragment "2" to fragment "1". It 
is argued that the difference between both boson transfer energies may serve as a criterion 
of how far we are from equilibrium. This difference is easily evaluated at each point of the 
multi-orbital energy surface and may speed up the search of the energy minimum. 

Explicit expressions have been determined within the frozen orbital approximation for the 
boson ionization potentials, boson affinities and boson transfer energies. These expressions 
shed light on the physical content of the quantities. The frozen orbital approximation is 
based on the assumption that adding or removing a single boson from a system with N 
bosons will not strongly change the orbitals of the other bosons and hence these orbitals 
can be assumed unchanged. We present numerical results which demonstrate that this 
approximation is valid for the condensates studied here. 

Finally, we would like to stress that the present findings are general and not at all re- 
stricted to the geometry of the double- well trap potential discussed here. We find similar 
results for other double- and multi-well trap potentials including potentials used in current 
experiments. 
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FIG. 1: (Color online) Results of one-orbital mean-field (GP) theory. Shown are the energies per 
particle Egp/N of the lowest energy solutions of the Gross-Pitaevskii equation as a function of 
A = \q{N —1). The bifurcation point observed at Acr = 0.837 is marked by a dash line. In the insets 
the GP orbitals corresponding to the upper (I) and lower (II) energy branches are depicted together 
with the double-well trap potential. The values of the external trap potential (for parameters, see 
text) have been scaled by 1/10. All energies are given in units of lo. All orbitals are dimensionless 
and plotted as functions of the dimensionless coordinate x. 
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FIG. 2: (Color online) Results of the two-orbital mean-field (MF(2)) theory for A = 0.8 < Xcr- 
Shown are energies per particle for condensates made of = 5, 10, 25, 10^, 10'^ and 10^ bosons as 
a function of the relative occupation number m/N of the orbital localized in the deeper well. The 
horizontal solid line labeled as "GP" shows the corresponding GP energy per particle. In the inset 
the energy per particle for the systems with = 13, 14 and 15 bosons are plotted. All energies 
are given in units of uj. 
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Occupation of the deeper- well orbital n^/N 

FIG. 3: (Color online) Results of the two-orbital mean-field (MF(2)) theory for A = 0.9 > Acr- 
Shown are energies per particle for condensates made of = 5, 10, 25, 10^, 10'^ and 10^ bosons as 
a function of the relative occupation number m/N of the orbital localized in the deeper well. The 
horizontal solid line labeled as "GP" shows the corresponding GP energy per particle. In the inset 
the energy per particle for the systems with N = 100000, 102000 and 105000 bosons are plotted. 
All energies are given in units of u). 
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FIG. 4: (Color online) Demonstration that BMF(2) provides the overall best mean field. Shown 
are the results obtained with the three-orbital mean-field (MF(3)) theory for = 1000 bosons and 
A = Ao(A^ — 1) = 0.9. Plotted are the energy per particle as a function of the relative occupation 
number ni/N for fixed values of the third relative occupation number n^/N. The curves shown 
correspond to nj,/N = 0.0, 0.0001, 0.001 and 0.005. All energies are given in units of oj. 
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FIG. 5: (Color online) Schematic diagram defining the boson ionization potentials (BIP) and boson 
affinities (BA). Total energies of condensates with N — 1,N and + 1 bosons obtained within the 
framework of the MF(2) theory are plotted in the left figure (A) as a function of the occupation 
number rii and in the right figure (B) as a function of the occupation number n2 = N — ni. 
The minima of the energy curves are denoted by ^~^Eq, ^Eq and ^^^Eq, and are indicated on 
the energy axis. The reference condensate for which the boson ionization potentials and boson 
affinities are defined is that with N particles. The optimal occupation numbers of this condensate 
are and n^. Vertical and adiabatic boson ionization potentials and boson affinities are shown 
(see text for more details). 
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FIG. 6: (Color online) The upper two panels show the chemical potentials (diagonal Lagrange 

multipliers) yUn and ^22 for condensates with = 25 and N = 1000 bosons as a function of 

the relative occupation number ni/N. For comparison, the horizontal solid line shows the value 

of the corresponding chemical potential /^gp obtained in the GP theory. The inset shows the 

same curves on a much (xlO^) enlarged scale. The lower two panels depict the derivative dE/dni 

(solid line) of the total energy with respect to ni and the difference between the boson transfer 

energies Q (points) for the same systems as in the upper panels. dE/drii is obtained by numerical 

differentiation of the MF(2) energy curves and Q is evaluated using Eg 1241 The vertical dashed 

lines mark the values of the optimal occupation numbers ni/N obtained at the minima of the 

corresponding energy curves. A = 0.9 is used throughout. The energy derivatives and all chemical 

potentials are given in units of w. 
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TABLE I: The adiabatic and vertical boson ionization potentials BIPa and BIPy{k) and boson 
affinities BA^ and BAyik) for condensates with N = 5,10,25,100 and 1000 bosons. Shown are 
quantities evaluated according to their definitions in Eas lllll'^ and I16I17I as differences of total 
energies and compared with those computed using the frozen orbital approximation in Easll4land 
1191 These approximate numbers are given in parenthesis. A = 0.9 is used throughout. 

Ionization Affinity 



N BIPa BIPv{l) BIPv{2) BAa BAv{l) BAv{2) 



5 -0.007508 0.025694 -0.006132 
( 0.027242) (-0.006248) 


-0.020167 -0.052911 -0.021267 
(-0.055309) (-0.021819) 


10 -0.010946 0.003733 -0.010390 
( 0.004460) (-0.010387) 


-0.016514 -0.031100 -0.017021 
(-0.031993) (-0.017142) 


25 -0.012613 -0.007123 -0.012412 
(-0.006844) (-0.012405) 


-0.014692 -0.020167 -0.014886 
(-0.020470) (-0.014910) 


100 -0.013361 -0.012032 -0.013314 
(-0.011964) (-0.013311) 


-0.013865 -0.015193 -0.013912 
(-0.015262) (-0.013916) 


1000 -0.013576 -0.013445 -0.013572 
(-0.013438) (-0.013571) 


-0.013626 -0.013758 -0.013631 
(-0.013764) (-0.013631) 



/iGP =0.013599 
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